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FRACTIONAL DIFFUSION IN GAUSSIAN NOISY ENVIRONMENT 


GUANNAN HU AND YAOZHONG HU 


Abstract. We study the fractional diffusion in a Gaussian noisy environment as described 
by the fractional order stochastic partial equations of the following form: Dfu{t, x) = 
Bu + u ■ , where is the fractional derivative of order a with respect to the time 

variable t, B is a. second order elliptic operator with respect to the space variable x G 
and a fractional Gaussian noise of Hurst parameter R = (Hi,-- - ,Hd). We obtain 
conditions satisfied by a and H so that the square integrable solution u exists uniquely . 


1. Introduction 

In recent years there have been a great amount of works on anomalous diffusions in the 
study of biophysics and so on (see for example [1], [4], [10], [11] to mention just a few). In 
mathematics, some of these anomalous diffusions (such as sub-diffusions) can be described by 
the so-called fractional order diffusion processes. As for the term “fractional order diffusion”, 
one has to distinguish two completely different types. One is the equation of the form 
dtu{t,x) = —(—A)"M(f, t), where f > 0, x G a E (0,2) is a positive number, ^ and 
A = Eti is the Laplacian. This equation is not associated with the anomalous diffusion. 
Instead, it is associated with the so-called stable process (or in general Levy process), which 
has jumps. Another equation is of the form D[°‘\{t,x) = Au{t,x), where is the 
Caputo fractional derivative with respect to t (see [9] for the study of various fractional 
derivatives). This equation is relevant to the anomalous diffusion which we mentioned and 
has been studied by a number of researchers, see for example [2] and the references therein. 
If one considers the anomalous diffusion in a random environment, then it is naturally 
led to the study of a fractional order stochastic partial differential equation of the form 
Uj“^M(t,x) = Bu{t,x) + u{t,x)W{t,x), where R is a second order differential operator, 
including the Laplacian as a special example, and IF is a noise. In this paper, we shall study 
this fractional order stochastic partial differential equation when W{t,x) = W^{x) is a time 
homogeneous fractional Gaussian noise of Hurst parameter H = {Hi, • • • , Hd). Mainly, we 
shall find a relation between a and H such that the solution to the above equation has a 
unique square integrable solution. 

If a is formally set to 1, then the above stochastic partial differential equation has been 
studied in [5]. So, our work can be considered as an extension of the work [5] to the case 
of fractional diffusion (in Gaussian noisy environment). Let us also mention that when we 
formally set a = 1, we recover one of the main results in [5]. Thus, our condition (2.10) 
given below is also optimal. 
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Here is the organization of the paper, Section 2 will describe the operator B and the noise 
and state the main resnlt of the paper. In onr proof we need to use the properties of 
the two fundamental solutions (Green’s functions) Z(t,x,^) and Y{t,x,^) associated with 
the equation D[°'^u{t^x) = Bu{t,x), which is represented by the Fox’s H-function . We shall 
recall some most relevant results on the H-function and the Green’s function Z{t,x,^) and 
Y{t,x,^) in Section 3. A number of preparatory lemmas are needed to prove main results 
and they are presented in Section 4. Finally, the last section is devoted to the proof of our 
main theorem. 


2. Main result 


Let 






dxndxj 




i=i 


be a uniformly elliptic second-order differential operator with bounded continuous real-valued 
coefficients. Let Uq be a given bounded continuous function (locally Holder continuous 
if d > 1). Let {W^{x) ,x G be a time homogeneous (time-independent) fractional 
Brownian held on some probability space (G,J^,P) (Like elsewhere in probability theory, 
we omit the dependence of W^{x) = W^{x,uj) on cj G G). Namely, the stochastic process 
{W^{x) ,x G is a (multi-parameter) Gaussian process with mean 0 and its covariance 
is given by 


( 2 . 1 ) 




i=l 


where Hi, • • ■ , Hd are some real numbers in the interval (0,1). Due to some technical diffi¬ 
culty, we assume that Hi > 1/2 for alH = 1, 2, • • • , d; the symbol E denotes the expectation 
on (D, B , P) and 

RHi{xi,yi) = ^ - \xi G M 


is the covariance function of a fractional Brownian motion of Hurst parameter H^. 

Throughout this paper we hx an arbitrary parameter a G (0,1) and a hnite time horizon 
T G (0,oo). We study the following stochastic partial differential equation of fractional 
order: 


( 2 . 2 ) 


D['^\{t,x) = Bu{t,x) +u{t,x) ■W^{x), tG(0,T], x G ; 

u{0,x) = Uo{x ), 


where 




7(1 — a) [dt 


d 

— / (t — r) °‘u{T,x)dT — t "m(0,x) 


is the Gaputo fractional derivative (see e.g. [9]) and W^{x) = (x) is the distri¬ 

butional derivative (generalized derivative) of , called fractional Brownian noise. 

Our objective is to obtain condition on a and H such that the above equation has a unique 
solution. But since is not differentiable or since W^{x) does not exist as an ordinary 
function, we have to describe under what sense a random held {^(t, x), t > 0, x G M'’*} is a 
solution to the above equation ( 2 . 2 ). 
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To motivate our definition of the solution, let us consider the following (deterministic) 
partial differential equation of fractional order with the term u{t, x) ■ W^{x) in ( 2 . 2 ) replaced 
hj fit, x): 


(2.3) 


D['^'^u{t,x) = Buit,x) + fit,x), tG(0,T], x G ; 
m ( 0 , x ) = uoix), 


where the function / is bounded and jointly continuous in (t, x) and locally Holder continuous 
in X. In [2] , it is proved that there are two Green’s functions 

{Z(t,x,0 , Yit,x,^) ,0 <t <T ,x,^ e M'’*} 


such that the solution to the Cauchy problem (2.3) is given by 

ds / Yit - s,x,y)fis,y)dy. 

jRd. 

In general, there is no explicit form for the two Green’s functions {Zit,x,^), H(f,x,^)}. 
However, their constructions and properties are known (see [2], [ 6 ], [7], and the references 
therein). We shall recall some needed results in the next section. 

From the classical solution expression (2.4), we expect that the solution uit,x) to (2.2) 
satisfies formally 


(2.4) 


uit,x)= / Zit,x,i)uoiCjdi + 


uit,x)=l Zit,x,^)uoi^)d^+ I ds I Yit-s,x,y)uis,y)W^iy)dy. 


The above formal integral Jq ds J^^Yit — s,x,y)uis,y)W^ iy)dy can be defined by Ito- 
Skorohod stochastic integral fo^it ~ s, x, i/)m(s, i/)ds W^idy) as given in [5]. 

Now, we can give the following dehnition. 


Definition 2.1. A random held |M(f, x), 0 < f < T, x G is called a mild solution to the 
equation (2.2) if 

(1) uit,x) is jointly measurable in t G [0,T] and x G M'’*; 

(2) V(t, x) G [0, T] X Jq Yit — s, X, y)uis, y)dsW^idy) is well dehned in 

(3) The following holds in 

(2.5) uit,x)=[ Zit,x,^)uoi^)d^+ [ [ Yit - s,x,y)uis,y)W^idy)ds. 

Jud. Jo jRd 

Let us return to the discussion of the two Green’s functions {Z(t,x,^), H(t,x,^)}. If 

d 

a = 1, namely, if in (2.3) is replaced by dt and B = A := then 

i=l 

(2.6) Zit, x,0 = Y it, X, 0 = (dTit)""*/^ exp 
In this case the stochastic partial differential equation of the form 

(2.7) ^ Auit, x) + u ■ W^ix), X G M'^, 
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was studied in [5]. The mild solution to the above equation (2.7) is proved to exist uniquely 
under conditions 

d 

(2.8) Hi >1/2, i = l,■■■,(! and Hj > d — 1. 

i=l 

The main result of this paper is to extend the above result in [5] to our equation (2.2). 


Theorem 2.2. Let the coefficients aij{x), hi{x) ,i,j = 1, • • • ,d, be bounded and continuous 
and let them be Holder continuous with exponent 7. Let aij{x) be uniformly elliptic. Namely, 
there is a constant oq G (0 , 00 ) such that 

d 

Y, > aoiep V e = (^1, ■ ■ ■ , ^d) e M'". 

*j=i 


Let Mo be a bounded continuous (and locally Holder continuous if d > 1). Assume 


(2.9) 
and 

( 2 . 10 ) 



if (i= 1,2,3,4 

if d > 5 


d 

YHi>d 

i=l 



a 


Then, the mild solution to (2.2) exists uniquely in 


Remark 2.3. (i) When a is formally set to 1, the condition (2.10) is the same as the 
condition (2.8) given in [5]. So, in some sense our condition (2.10) is optimal. 

(ii) Since Hi < 1 for alH = 1,2, ■ ■ ■ , d the condition is possible only when a > 1/2. 


3. Green’s functions Z and Y 


3.1. Fox’s H- function. We shall use H-function to express the Green’s functions Z and Y 
in Dehnition 2.1. In this subsection we recall some results about the H-function s and the 
two Green’s functions. We shall follow the presentation in [8] (see also [2] and references 
therein). 


Definition 3.1. Let m, n,p, q be integers such that 0<m<q',0<n<p. Let Oj, 6* G C be 
complex numbers and let aj, /3j be positive numbers, i = 1, 2, • • • ,p; j = 1, 2, • • • ,q. Let the 
set of poles of the gamma functions r{bj + fdjs) doesn’t intersect with that of the gamma 
functions r(l — a* — Ojs), namely. 


for all i 



~bj — 1 — ai + k 

, p and j = 1,2, ■ ■ ■ ,q. The H-function 



Trmnt — TTvan 

Ppq [^) = Ppq 



(ai,ai) ■■ 


z 

(&i,A) •• 

■ (bq.Pg) _ 


is dehned by the following integral 


(3.1) 


1 f n7im+fts)nr..r(i-oi-M _ 

a'r* Jl nL.+. r(s + ais) r(i - b, - i3,s) ‘ “ 


2 ( G C, 
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where an empty product in (3.1) means 1 and L in (3.1) is the infinite contour which separates 
all the points bji to the left and all the points Qik to the right of L. Moreover, L has one of 
the following forms: 

(1) L = L_ao is a left loop situated in a horizontal strip starting at point —oo + icpi and 

terminating at point — cx) + for some — oo < < 02 < oo 

(2) L = L_|_oo is a right loop situated in a horizontal strip starting at point +oo + i0i and 
terminating at point oo + i 02 for some — oo < 0i < 02 < oo 

(3) L = Li^oo is a contour starting at point 7 — ioo and terminating at point 7 + ioo for 
some 7 G (— 00 , 00 ) 

To illustrate L we give the following graphs. 


/ 

\ 



X 



• . 0 


- 7 - 

Case 1 




Case 2 


Case 3 


The integral (3.1) exists when Yl'j=i l^j ~ Y7i=i ^ 0 (see [ 8 ], Theorem 1.1). 

Example 3.2. To compare with the classical case a = 1, we consider the case m = 2, n = 0, 
p = 1, g = 2, oi = «! = 62 = /Si = /d 2 = 1 and hi = Let L = L_^. Then, we have 


H 


( 1 . 1 ) 

(f.i). (1.1) 


(3.2) 


1 

2'Ki 


r(i + ")r(i + T-* 


r(i + s) 


d 


V—\ , d , ,(i 

= 2 ^ hm (s +- +n)r(-+ s)^ 


i ;=0 

00 


lim 


r(n + | + s + l) 


r\ ^^ V o 
V = \) ' 2 


-(i+v) (s + I + n- l)---(s + f) 

1 
v\ 


yd/2 


= z / exp(- 2 :). 


3.2. Green’s functions Z and Y when B has constant coefficients. In this subsection 
let us consider Z and Y when the operator B in (2.2) has the following form 




02 


Li=i 


dxjdxi ’ 
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where the matrix A = (aij) is positive definite. In this case, Z and V (we call them Zq and 
Yq to distinguish with the general coefficient case) are given as follows. 

-d/2 r AL 1 


Zo(t,x) = 


TT 


(det Ay/'^ 




/■) 


yc 2 X j 




xHll 


d 


i,j=l 


(l,a) 

( 1 , 1 ) 


where = A ^ and 


Yo{t,x) = 


TT 


-d /2 


r d 


(det Ay/"^ L 


E^' 


U) 


Ob'iOb j 


1 -d /2 

t 


OL—1 


x-ff?? 




E-^' 


u) 


ObiOb j 




(a, a) 

(i,i), (1,1) 


It is easy to see that for the constant coefficient case, the both of the Green’s functions are 
homogeneous in time and space. Namely, 

Zoit, X, 0 = zyt, x-i), Yo{t, X, 0 = Yo{t, X - 0 • 

In particular, when a = 1, it is easy to see from the above expression and the explicit form 
(3.2) of ^^ 12 ( 2 ^) that 

E?,=i - ^i){xj - Q 


Zo(t,x,0 = Yo{t,x,^) = (dvr) '^/Met(Al) 


exp 


4t 


which reduces to (2.6) when A = J is the identity matrix. 

With the above expression for Zq and Yq and the properties of the fl-function , one can 
obtain the following estimates. 


Proposition 3.3. Denote 


(3.3) 


p(t,x) = exp ( — at 2 -a|a;| 2 -a^^ t>0. 


X G 


where a G (0, cx)) is a constant whose exact value is irrelevant in the paper. Then, we have 
the following estimates: 


(3.4) 


Ct 2 py,x) when d = 1 

|^o(^, a;)| < ■{ log ^1 + l]p(t, x) when d = 2 

Ct~'^\xy~'^p{t,x) when d>3, 


where for instance, \Zo{t,x)\ < Ct~'^p{t,x) means that there are positive constant C and 
positive constant a such that the above inequality holds. In what follows the positive con¬ 
stants C and a are generic, which may be different in different appearances. 

Proof. Denote R = x^/t". From [2], Proposition 1, it follows that when i? < 1, we have 

{ Ct-1 when d = 1 

Gt-“[|log^|+ 1] when d = 2 
Ct-yx^-'^ when d>3. 
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Since when R < 1, p{t,x) is bounded from below. This proves the inequality (3.4) when 
R<1. 

When R> 1, then by [2], Proposition 1 we have |Zo(t, x)| < Ct~^p(t,x). It is clear that 
this implies the inequality (3.4) when d = 1 and d = 2. Now, we assume that d > 3. We 
have 

\Zo{t,x)\ < cr’^p{t,x) <cr°‘x^~‘^ pit,x) 

< Ct-‘^x^-^p{t,x), 

( 2\f-l 

where we used the fact that ( ^ j p{t, x) < p{t, x) for a different a in the later p{t, x). □ 

Similarly, we can use [2], Proposition 2 (for d = 1 case) and [2], Section 4.2 (for d > 2 
case) to obtain the following estimates for YQ(t,x). 


Proposition 3.4. We follow the same notation p{t,x) as dehned by (3.3). We have 

(1) When d = 1, we have the following estimates: 

(3.5) |>o(f,a;)| < 

( 2 ) 


(3.6) 


Cf 2 ^p[t,x) when t “|a:f > 1 


when < 1. 

When d > 2, we have the following estimates: 


\Yo{t,x)\ < < 


Ct ^p{t,x) 

when 

d = 2 

Ct~^~^p{t, x) 

when 

d = 3 

[| log ^ 1 -f l]p(t, x) 

when 

d = 4 

Ct~°'~^\x\^~'^p{t, x) 

when 

d > 5, 


where for instance, |yo(f)3;)| < Ct~^p{t,x) means that there are positive constant C 
and positive constant a such that the above inequality holds. In what follows the 
positive constants C and a are generic, which may be different in different appear¬ 
ances. 


3.3. Green’s functions Z and Y in general coefficient case. If the coefficients of B 
are not constant, then the Green’s functions Z and Y are more complicated and may be 
obtained by a method similar to the Levi parametrix for the parabolic equations. 

Denote 


d 




M{t,x,^) = ^[aij{x) - aij{^)]-^^^Zo{t,x - ^) 

i,j=l ® ^ 


d 


d 


+ ^bi{x) — Zo{t,x - 0 + c{x)Zo{t,x- 0 


2 = 1 




i,j=l ® ^ 


d 


d 


+ ^bi{x)—Yo{t,x - + c{x)Yo{t,x - ^) 


2 = 1 


K{t,x,^) 
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Let Q{s,y,^) and d)(s,|/, be defined by 

Q{t,x,0 = M{t,x,0 + [ ds [ K{t - s,x,y)Q{s,y,^)dy, 

Jo 

= K{t,x,0 + [ ds f K{t - s,x,y)^{s,y,i)dy 
Jo Jr'^ 

Proposition 3.5. Let the coefficients aij{x) and bi{x) satisfy the conditions in Theorem 2.2. 
Recall that 7 is the Holder exponent of the coefficients with respect to the spatial variable 
X. Then, the Green’s fnnctions {Z{t,x,^),Y{t,x,^)} have the following form: 

Z{t, X, 0 = Zo{t, x-^) + Vz{t, X, 0; 


(3.7) Yit, X, 0 = Yoit, X-0 + Vvit, X, 0 , 

where ^ 

Vz{t,x,^)= ds Yo{t - s,x,y)Q{s,y,^)dy; 

Jo JR‘i 

yY{t,x,^)= ds Yo{t - s,x,y)^{s,y,^)dy. 

Jo Jr<^ 

Moreover, the fnnction Vzit,x,^),VYit,x,^) satisfy the following estimates. 


(3.8) 


\Vz{t,x,^)\ < 


'Ct^ 2p(f, X — , 

Ct^-^p{t,x - , 

Ct^ ^\x — X — ^), 

^(7t(7-70)f-«|a; _ ^|-2+7-270p(^^3. _ ^ 


and 


(3.9) 

\VY{t,X,^)\ < 


' Ct°‘ ^^2p(t,X — , 

Ct^-^p{t, X - , 

Gf(70+7)f-l|a; — X — , 


when d = 1 ; 
when d = 2 ; 
when d = 3 or d > 5 ; 
when d = 4 

when d = 1 ; 
when d = 2 ; 
when d = 3 or d > 5 ; 
when d = 4 


Here 70 is any nnmber snch that 0 < 70 < 7 and in the case d > 3, the constant C depends 
on 7 o. 


4. Auxiliary lemmas 

To prove onr main theorem, we need to dominate certain mnltiple integral involving 
Y(t,x,^) and Z(t,x,^). Since both Y(t,x,^) and Z{t,x,Q are complicated, we shall first 
bonnded them by p(f, X—^) from the estimations of |To(^,^) I and |Hy(f,x,^)|. More precisely, 
we have the following bonnds for Y(t,x,^). 

Lemma 4.1. Let x G G (0,T]. Then 

Ct~^~^‘5p[t,x — d=l; 

Ct~^p{t,x — d = 2 ; 

Ct~y~Yo)^~^\x — X — d = 4; 

( 7 ^-( 7 - 7 o)f-i|^ _ ^| 2 -d+( 7 - 7 o)/ 2 p^^^ ^ d = 3 ord> 5 . 


( 4 . 1 ) \Y{t,x,0\< 
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Proof. We shall prove the lemma case by case. First, when d = 1, hj Proposition 3.4, we 
have 




If t "|a; — < 1 , then 


Ct2 — t "lx —> 1 ; 

< 1 . 


,_i+ 5 | p{x,t) 


\Yo{t, X - f)\ < Ct < Ct^ p{t,x-f). 

Therefore 

|>"(i,x,OI < |>o(i,x -01 + 

< X - 0 + C'f-^+tp(t, X - 0 

< X — 0 • 

Now we consider the case d = 2. From the following inequalities: 

|W^(i,x,OI < CP'^-^p{t,x-f); 
l^o(^,a:-OI < Ct-^p{t,x - f) 

we have easily 

|h'(t,x,OI < |>o(t,x-OI + l'^y(^,a;,OI < Ct-^p{t,x - f). 

We re going to prove the lemma when d = 3. From Proposition 3.4 we have 
|Fo(t,x-OI < Ct~‘^~^p{t,x - f) 


— (^^-(7-7o)f-l|2; _ ^|-1+(T'-7 o)/2 


X — f 


t 2 


1-(7-7o)/2 


p(t,x -0 


< Ct — ^1 ^ + P X — 0 • 

Combining this inequality with Proposition 3.5 we obtain 

|r(t, X, 01 < X - 0 • 

We turn to consider the case d = 4. Proposition 3.4 yields that for any d > 0 the following 
holds true: 


|Fo(Ca:-OI < Ct 


— OL —1 


|x-0 


2 \ 6 


+ 








= ct 


— OL—l 


|x- 0 ^ 


e r 




2\ 26» 


+ 1 


p(t,x-0 


If ^0^ ^ 




+ 1 


p(t,x -0 < 2 




2 \ 2 e 


As a consequence, we have 

|ho(Cx-OI 




p(t,x-0 < C'p(t,x-0 


p(t,x- 0 - 
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If < 1 , then the above inequality is obviously true. Now, we can choose 9 > 0 , such 

that —29 > (—2 + 7 — 270). Thus, we have 

|Fo(t, 2 :-OI = |a;-^|- 2 +T'- 2 T'o 

_ f I \ - 26 »-(- 2 + 7 - 270 ) 

p{t,x-0 

< . 1 ^; _ X — . 

Combining the above inequality with Proposition 3.5 we have 

_|_(7^(70+7)f-l|3™ _ X — ^) 

since -(7 - 27o)f - 1 < (70 + 7)f - 1. 

Finally we consider the case d > 5. From the estimates |yo(t,a^ ~ 01 — Ct~°'~^\x — 
^|4-(ip(^, X — ^) we obtain 



|bo(t,2:-0l < 


x-i 


t2 


2-(7-7o)/2 


p(t,x-0 


< i-(7-70)f-l 1^ _ ^|2-d+(7-70)/2p('^^ X - 0 • 


Therefore, we have 

\Y (f, X, 01 < X — 0 • 

The proposition is then proved. □ 


The bound (4.1) will greatly help to simplify our estimation of the multiple integrals that 
we are going to encounter. However, when the dimension d is greater than or equal to 2, the 
multiple integrals are still complicated to estimate and our main technique is to reduce the 
computation to one dimensional. This means we shall further bound the right hand side of 
the inequality (4.1) by product of functions of one variable. Before doing so, we denote the 
exponents of t and |x — in (4.1) by Q and Hd- Namely, we denote 


(4.2) 


and 

(4.3) 


Crf 


'-1 + f, d 

— 1 , d 

-(7-27o)f-l, d 

,-(7-70)7-1, d 


1 ; 

2 ; 

4; 

3 or d > 5 . 


[ 0 , d = l, 2 ; 

Kd= 1-2 +'y- 270 , d = 4; 

1^2 — d + (7 — 7o)/2, d = 3 or d > 5 . 


From now on we shall exclusively use p{t,x) = exp ( — (jt“ 2 ^|x| 2 ^) to denote a function 
of one variable. However, the constant a may be different in different appearances of p{t, x) 
(for notational simplicity, we omit the explicit dependence on a of p{t,x)). 

With these notation Lemma 4.1 yields 
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Lemma 4.2. The following bound holds true for the Green’s function Y: 

d 

(4.4) \Y{t,x,0\<Cl[&/^\xi-fr‘^/^p{t,x,-f,). 

Proof. It is easy to see that 


2=1 


\x\ = 


d \ V2 
„2 


E xj 1 > max \xi\ > TT \xi\d 

I i<i<d J. J.' ' 


, 2=1 


2=1 


Thus for any positive number a > 0, |x| " < ntikr^- 
On the other hand, 


2 

x\ 2 -« 



*- 2=1 


max \Xi 

l<i<d 


1 

2-ck 


max 

i<2<d d ^ ^ 

2 = 1 


Combining the above with (4.1) yields (4.4) since the exponents in \x — in (4.1) are 
negative. □ 


Lemma 4.3. Let —1 < /? < 0, a; G M. Then, there is a constant C, dependent on a, a and 
(3, but independent of f and s such that 


sup / \x\^p{s, X — f)dx < Cs ^. 
CeM Jr 


Proof. Making the substitution x = ys 2 we obtain 


\xfp{s,x — f)dx = s 2~^2 


Oi/3 I a 

< S 2 "^2 


• exp — a 


S2 


dy 


\y\^dy + 


<1 


/ exp ( — a 

y —ET 

Jr V 

S2 


2 

2-a. 


dy 


< CS 


since the two integrals inside the parenthesis are finite (and independent of s and f). 
The following is a slight extension of the above lemma. 


□ 


Lemma 4.4. There is a constant C, dependent on a, a and (3, but independent of f and s 
such that 
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Proof. We shall follow the same idea as in the proof of Lemma 4.3. Making the substitution 
X = ys 2 we obtain 


\xf I log |x| I p(s, X — f)dx 

[ |i/|^[|log|i/|| + |logs|]-exp ( - a 


y- 


52 


2 

2-a 


dy 


< C'sf + 2(i +|iogs|)( / |i/|^|log|2/||d|/+ / exp ( -a 

\ J |t/|<e J M 

< (1 + I log s|). 


y 


S2 


dy 


This proves the lemma. 


□ 


Lemma 4.5. Let 9i and 6*2 satisfy —1 < 6*1 < 0, —1 < 02 < 0. Then for any pi, T 2 G M, pi 7 ^ 

r2, 

(1) If 9i + 02 = —1; then 


\pi - - Pif^p{s2 - SI,P2 - Pi)dpi <C + C\ log(p2 - Ti) 


(2) 7/01 + 02 < —1, then 


\pi - 'ri|^^|p2 - Pif''p{s2 - si,p2 - pi)dpi < C\p2 - 


Proof. Without loss of generality, suppose Ti < P 2 - We divide the integral domain into four 
intervals. 


|pi - 'rir"|p2 - Pi|%(s2 - si,p2 - pi)dpi 


2 


' —00 


|Pl - ^l|^'|P2 - Pir"p(s2 - Si,P2 - pl)dpl 


| 6»2 


+ / |pi - ri|®l|p2 - Pl|'^^p(s2 - S1,P2 - Pi)dpi 

/ 3ti-p2 

-2- 

3p2-n 

+ / |Pl - Ui|®l|p2 - Pl|'^"p(s2 - Si,p2 - Pl)dpi 

J ^1+P 2 

poo 

+ / |pi -'ri|®i|p2 - Pi|%(s2 - S1,P2 - Pi)dpi 

/ 3P2-Ti 
-2- 


—: Ji + J2 + Is + h 
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Let US consider I 2 first. When pi G 
p{s 2 — si, p 2 — pi) < 1 , we have the 


3ti — P 2 Ti + P 2 


2 ’ 2 
ollowing estimate for l 2 '- 


, we have \p 2 — pi| > • Noticing 


h < 


< 


02 pn±£i. 

P 2 - n \ r 2 


P2 -Ti 


O 2 


' 3ti-P2 
2 

f 2 


|pi - Til'^Mpi 
(Pi - 'ri)'^Mpi + 


= C* P2 - Ti 


T 1 

1 +01+02 


f 3 ti-P 2 
2 


(ti - pi)^Mpi 


With the same argument, we have 


h<C P2-T1 


1+01+02 


Now, we study Ji. The term 14 can be analyzed in a similar way. Since pi < < 

Ti < P 2 , we have 


, 3ti-p2 


/l< 


+ - - -Si, p 2 - pl)dpi 


To estimate the above integral, we divide our estimation into three cases. 
Case i): 6*1 + 6*2 < —1. 

In this case, we bound p(s 2 — -Si, p 2 — pi) by 1. Thus, we have 


Ii< 


, 3ti-P2 
2 


(ti - pif^^^^dp, = 


1 + 61 + 62 


P2 - Ti 


1+01+02 


Case a): 61 + 62 = —1, > L 

In this case, we have < Ti — 1. Thus, we have 


h < 
< 

< 


r-Tl-l 


(ri-pi) p(-S 2 --Si,P 2 - Pi)dpi 


' —00 

pTi-1 


’ —00 
poo 


P{S2 - Si,p2 - Pl)dpi 
P{S2 - Si,p2 - Pl)dpi 


which is bounded when si and S 2 are in a bounded domain. 
Case Hi): 61 + 62 = —1, < 1- 
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In this case, we divide the integral into two intervals as follows. 


h = 


■ 3ti-P2 
2 


('Tl - - Si, p2 - Pl)dpi 


' —oo 

PTl—1 


< / {Ti- pl)~^p{s2- Si,p2- pl)dpi+ in- Pl)~^p{s2- Si,p2- Pl)dpi 

J — OO j Ti — 1 

<C+ {Ti-pi)~^dpi 

J Ti — 1 

< C' + C'|ln(/i 2 -ri)|. 

Similar argument works for 14 . Combining the estimates for Ik,k = 1,2, 3,4 yields the 
lemma. □ 

Lemma 4.6. Let 9i and 62 satisfy —1 < 6^1 < 0, —1 < ^2 < 0 and 9i + 202 > —2. Let 
0 < Ti < r 2 < T and 0 < Si < S 2 < T. Then for any pi, T 2 G M, pi 7 ^ T 2 , we have 


\pi - Til 1|P2 - Pl| ^17-2 - Til ^p{s2 - S1,P2 “ Pl)p{r2 “ ri,T2 - Ti)dpidTi 

2 

C(S2 - Si)''^'^(r2 - n)^, 01 + 02 > -1; 

(4.5) < <; C(r2-ri)^'^^^, 0i + 02 <-1; 

C(r 2 - [1 + I log(r 2 - ri)|] , 0i + 02 = -1. 

Proof. First, we write 


I : = 


\pi - nf"\P 2 - Plf^\r 2 - - Si,p 2 - Pl)p{r 2 - ri,T 2 - Ti)dpidTi 


1^2 


(4.6) = / /(ri,p2,si,S2,0i,02)|r2 - - ri,r 2 - ri)dri, 


where 


/(ri,p2,Si,S2,0i,02) = / |pi-ri|®i|p2 - Pir"p(s2 - Si,p2-Pi)dpi. 


I® 2 „ 


We divide the situation into three cases. 

Case 1 ): 01 + 02 > —1. 

In this case we apply the Holder’s inequality to obtain 


/('ri,P2,Si, 52,01, 02) < i / |pi - 'ri|®i+^’'p(s2 - S1,P2 - Pi)dpi 


G 

6 ' i +^2 


|P2 - Plf"^^^p{s2 - SI,P2 - Pl)dpi 


^2 

^1+02 


(4.7) 


^(61+82) I a 

< C (s2 — Si) 2 2 
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where the last inequality follows from Lemma 4.3. Substituting the above estimate (4.7) into 
(4.6), we have 


/ = 


/(T1,P2, Si,S2,9i,92)\t2 - ri|%(r2 - ri,T2 - Ti)dTi 

|t 2 - ri|%(r 2 - ri, Ts - Ti)dTi. 


a( 6 -\ + 0 o) , a 

< C{S2-Sl) 2 +2 


Using Lemma 4.3 again we have, 

ce(0-i -j-Oo) , a Oi9o , a 

I < C{s2 — Si) 2 +2 (r2 — ri)^'''2 . 

Case a ): + 6*2 < —1. 

In this case, from Lemma 4.5, part (ii) it follows 

/(ti,P2,si,S2,6'i,6'2) < C\p2 - ri|®i+®2+i _ 

Hence, we have 

I < C f \p2-Tif^^^^^^\T2-Tif^p{r2-ri,T2-Ti)dTi. 

Jr 

Now, since from the condition of the lemma, 6^1 + 202 + 1 > —1, we can use Hblder’s inequality 
such as in the inequality (4.7) in the case (i), to obtain 

I < C [ \p2-Ti\^^^^^^'^\T2-Tif^p{r2-ri,T2-Ti)dTi 

Jr 

q;( 0 -| +26o) I 

< C(r 2 — ri) 2 . 

Case iii): 9i + 92 = —1. 

In this case, we hrst use Lemma 4.5, part (i) to obtain 

/(U,P 2 ,Sl, 52,01,02) < C [1 + |log|p 2 - Till] . 

Thus, using Lemma 4.4, we have 

I < C f{l + \\og\p 2 -Ti\\}\T 2 -nf^p{r 2 -ri,T 2 -Ti)dri 


< C(r 2 -ri) 2 [1 + |iog |r 2 - ri| 


The lemma is then proved. 


□ 


Corollary 4.7. Let 0i and 02 satisfy —1 < 0i < 0, —1 < 02 < 0 and 0i + 202 > —2. Let 
0 < n < r 2 < T and 0 < si < S 2 < T. Then for any pi, r 2 G M, pi 7 ^ T 2 , we have 


< 


|P1 - Til ^\p2 - Pl\ ^|t 2 - Til ^P{S2 - Si,p2 “ Pi)p(t2 “ ri,T2 - Ti)dpidTi 

^ a(ex+292+2) ^ ^ 0(91+292+2) 

C(S 2 - Si) 4 (r 2 - ri) 4 ; 

C(S 2 - Si)^^^^(r 2 - [1 + I log(r 2 -ri)| + | log(s 2 - Si)|] ; 0^ + 02 = -1. 


01 + 02 7^ —1 


(4.8) 


Proof. Consider first the case 0i + 02 < —1. Denote the integral on the left hand side of (4.8) 
by I. Then the inequality (4.8) implies 

cx.(0-\ + 200 ) I 

/ < C{r 2 -ri) —^—+“. 
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In the same way we have 


Q!(6>i+2g2) 


+a 


I < C(S2 - Si)" 

Now we use the fact that if three numbers satisfying a < b and a < c, then a = < 

51 / 2 ^ 1 / 2 _ 

I < C{r2 - 

which simplihes to (4.8). The exactly the same argument applied to the case 6 *i + 6*2 = — 1 
and the case 6*1 + 6*2 > —1. Thus, the inequality (4.7) implies (4.8). □ 


Lemma 4.8. Let pi, ■ ■ ■ > 0. Then for any T > 0, 


r"nLir(pO 


(4.9) / (s„-s„_i)'’" ‘■■■(S 2 -S 1 )” 'sf ‘<is = 

Jo<si<-<Sn<T T{pi + ■ ■ ■ + Pn + 1) 

Proof. This is well-known. For example it is straightforward consequence of formula 4.634 
of [3] with some obvious transformations. □ 

Lemma 4.9. Assume that Uq is bounded. Then 


sup / Z{t,x,^)uo{0d^ <C. 
xsM Jm<i 

Proof. We use Z(t, x, = Zo(t, x — ^) + Vz(t, x, ^). Since ug is bounded, 


Zo(t,x,^)uo( 0 d^ 


< C / jZo(t,x,Oid^ 


which is bounded by the estimates in (3.4) and a substitution ^ = x Pt^y. In fact, we have, 
for example, when d > 3, 

f \Zo{t, X — ^)\d^ < C f t~°‘t^~^t^\y\‘^~'^exp{—a\y\^^}dy<Ct^~°‘<C. 

JR'* JRd 

Similarly, using the estimation for Vz(t,x,^) given in Proposition 3.5 we can bound 
/iRd \ Vz(t,x,^)\d^ by a constant. In fact, for example, when d = 3, we have 

' - / 3a (t —1 2 'ya 

/ t^t 2 exp{—cT|?/|^}(ii/< < C. 


/ \Vz{t,x,0\d^ < Ct--^ 

JR'* jR‘i 

The other dimension case can be dealt with the same way. 


□ 


5. Proof of the main theorem 2.2. 

Change t to s and x to y and the equation (2.5) for mild solution becomes 

u{s,y)= f Z{s,y,0uo{0d^+ [ [ Y{s - r,y, z)u{r, z)W^{dz)dr . 

jRd- Jo JRd 

Substituting the above into (2.5), we have 

u{t,x) = [ Z{t,x,0uo{^)d^+ [ [ Y{t - s,x,y)Z{s,y,Ouo{Od^W^{dy)ds 

jRd, Jq J]R2d 

+ f f f Y{t — s,x,y)Y{s — r,y,z)u{r,z)W^{dz)drW^{dy)ds . 

Jo Jo iR2'i 
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We continue to iterate this procedure to obtain 

OO 

(5.1) u{t,x) = , 

n=0 

where satishes the following recursive relation: 

= [ Z{t,x, 0 uo{ 0 d^ 

rHt 


'h„+i(t,x) = / / Y{t-s,x,y)'^n{,s,y)W {dy)ds, n = 0,1,2, 

Jo 

To write down the explicit expression for the expansion (5.1), we denote 

(5.2) fn{t,x-,Xi, - ■ ■ ,Xn) = / / Y{t-Sn,X,Xn)---Y{s 2 -Si,X 2 ,Xi)Z{si,Xi,^)UQ{^)d^ds, 

Jt„ Jr'* 

where 

Tn = 0 < Si < S 2 < ■ ■ ■ < Sn < t and ds = dsids 2 ■ ■ ■ dsn ■ 

With these notations, we see from the above iteration procedure that 

T„(t,x) = Inifnit,x)) 

= [ fn{t,X-,Xi,- ■ ■ ,Xn)W^{dXi)W^{dX2) ■ ■ -W^idXn) 


(5.3) 


fn{t, x; Ti, • • • , Xn)W^ {dXi)W^ {dX2) ■ ■ ■ hh (dXn) ■ 


where In denotes the multiple Ito type integral with respect to W (x) (see [5]) and fn{t, x;xi, - ■ ■ 
is the symmetrization of fnit, x;xi, - ■ ■ , Xn) with respect to xi, • • • , x„: 

fniji ^ 1 ' ' ' ) ^n) T ^ ^ fn(ti X, Xj^ , • • • , Xj ) , 

nl 

,Inherin') 

where a(n) denotes the set of permutations of (1,2, • • • ,n). 

The expansion (5.1) with the explicit expression (5.3) for \I/„ is called the chaos expansion 
of the solution. 

If the equation (2.2) has a square integrable solution, then it has a chaos expansion ac¬ 
cording to a general theorem of Ito. From the above iteration procedure, it is easy to see that 
this chaos expansion of the solution is given uniquely by (5.1)-(5.3). This is the uniqueness. 

If we can show that the series (5.1) is convergent in J^, P), then it is easy to verify 

that u(t,x) dehned by (5.1)-(5.3) satishes the equation (2.5). Thus, the existence of the 
solution to (2.2) is solved and the explicit form of the solution is also given (by (5.1)-(5.3)). 
We refer to [5] for more detail. 

Thus, our remaining task is to prove that the series dehned by (5.1) is convergent in 
P). To this end, we need to use the lemmas that we just proved. 
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Let 

now u(t,x) be defined by (5.1)-(5.3) 

Then 

we have 



E,[u{t,x)‘^] = 

00 

n =0 

In{fn{t,x)) 

2 


= 

00 

^n\{fnJn)H 

n =0 


(5.4) 

< 

00 

n =0 


where 





(5.5) 

P n 

J^2nd 


)h(ni,--- W, 

T^)duidvidu 2 dv 2 


and the last inequality follows from Holder inequality. Here and in the remaining part of the 
paper, we use the following notations: 


Xli ‘ ‘ 5 dUi dUi^ * * * dUid , t 1 , 2 ,***, 77 ', 

d d 

Vi) = JJ (p//. (Uij, Vij) = JJ Hj {2Hj - 1) \Uij - Uy . 
i=i i=i 

We use the idea in [5] to estimate each term Qn(t,x) = n\{f,f)H in the series (5.4). By 
the defining formula (5.2) for /„ we have 


Qn{t,x) = n\ / 

Jt2 7R2nd+2 

■ / Z{si, ^ 1 , ^o)uo{^o)d^o • y (f - r„, x,rin)---Y (ra - n, 773, hi) 

jRd 

■ / Z{ri,r]i,r]o)uo{r]o)dr]od^dT]dsdr. 

Jr'^ 

Application of Lemma 4.9 to the above integral yields 

Qn{t,x) < Cn\ I / Y\ - Vi)Y{t - Sn,X,^n) ■■ ■Y{s 2 - Si,^2,^1) 

Jt 2 J^2nd 

■Y {t - Tn, x,r]n)---Y (ra - n, r/a, r]i)d^dr]dsdv. 


Using Lemma 4.2 to the above integral, we have 
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where 


Qi^n(t,Xi, s,r) = 


/ \Tl^HXpk-rk)>\t 

U=i J 


I 


Sn\ \Xi - Pn\ P{t - Sn,Xi - Pn) 


Cd 

• ■ |S2 - Si|^|p2 - Pl\~^P{s2 - Sl,P2 - Pi) 

\t - rn\^\Xi - Tn\^p{t - rn,Xi - T^) • • • |r 2 - n 
, Id 




■|t 2 - Til d p{r 2 - ri,T 2 - Ti)dpdT. 

Here we use the notation p^ = ^ki and = pki) k = 1, • • • ,n. The quantity can be 
written as 


Cd , 


0 i,„(t,Xi,s,r) = \t-Sn\d\t-r, 


I Cd 


I I 


I 

d 


/ \\\p>HApk-rk)\ 

U=i J 


S2 - si| d \r2 - ri| 

Kd 

- Pn|“P(i - Sn,Xi - p„) 


(5.7) 


Kd '^d 

■\Xi - Tnl^pit - rn,Xi - T^) ■ ■ ■ \p2 - Pl\^ p{s2 “ Si, P2 - Pi) 

Kd 

• 17-2 - ri|^p(r2 - ri,r2 - Ti)dpdT . 


From the dehnition (4.3) of Kd we see easily ^ > — 1. We assume 
(5.8) 


2 Hi + ‘^> 0 . 

d 


Under the above condition we can apply the Corollary 4.7 with 6 i = 2Hi — 2 > —1, 62 = 
^ > —1 to the integration dpidri in the expression (5.7) (Condition (5.8) implies that 
9i + 26*2 > —2). Then, when 6*1 + 6*2 7 ^ —1, we have 


<d 


Qi^n{t,Xi,s,r) < C\t - Sn\ d \t - r. 


I ^d 


■ S2 - Si 


, Cd I Hjd+I^d , 

2d 


I I — I 

|S3 - S2I |r3 - r2 

Cd I gjd+Krf 

\T2 — 'k'l\ d ^ 2d 


^d 

d 


/ \TlPHXpk-rk)\ 

U=2 J 


Kd 

\Xi - pn\^p{t - Sn, Xi - pn) 


Kd Kd 

■\Xi - Tn\^p{t - rn,Xi - T^) ' ' ' |p3 “ P2|^P(S3 “ S 2 , P 3 “ P 2 ) 

Kd 

■ ■■ 17-3 - nl^p^rs - r 2 , n - T 2 )dpn ■ ■ ■ dp 2 dTn ■■■dT 2 . 
Repeatedly applying this argument, we obtain 

n 

(5.9) Qi ^ n { t , Xi , s , r ) < C” JJ|4+i-4|^‘|sfc+i-Sfcl^S 

k=l 


where we recall the convention that tn+i = t and s„+i = s and where 


n Cd Hid -\- Kd 

"i = — H-- a 


d 


2d 
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Substituting the above estimate of Qi^n into the expression for 0„, we have 

en{t,x) < C'’^ 


= C" 


/ It 

]^(sfe+i - SkY{rk+i - VkYdsdr 

n k = l 

„ n ' ^ 

I ('^fc+i Sk') ds , 


k=l 


where 


— Cd + 


+ with = 

9 9 ' ' ^ * 


Z =1 


1=1 


Now, we apply Lemma 4.8 to obtain 

Qn{t,x) < 


< 


r(£ + i) 
_r(n(£ + l))_ 

cn 


r(2n(£+l)) ■ 

This estimate combined with (5.4) proves that if 

(5.10) 2(£ + l)>l, 

then YIYLo 0n(t, x) is bounded which implies that the series (5.1) is convergent in L^(fl, J^, P). 

From the explicit expressions of (d and Kd we see by analyzing the condition (5.10) for the 
cases d=l, d = 2, d = 4: and d = 3 or d > 5 separately. We see that the condition (5.10) is 
equivalent to 


(5.11) 


^ ^ Hi > d — 2 + 


i=l 


a 


When 6 i + 02 = —1, Corollary 4.7 implies that for any £ > 0, 


/ \pl - n|^|P 2 - Pi| "|t -2 - riY'p{s 2 - S 1 ,P 2 - Pi)p{r 2 - ruT 2 - Ti)dpidTi 

d]R2 

, a;(6>2 + l+£) , , c,(e2 + l++s) 

< C{S 2 - Si) 4 (r 2 - n) 4 

Now we can follow the above same argument to obtain that if 

(5.12) 2(£ + l)>l, 

where t = then 0n(t, x) is bounded. In the same way as in the case 0\P 02^ “I; 

we can show that the condition (5.11) implies (5.12). 

Now we consider the condition (5.8). From the dehnition (4.3) of Kd, we see that when 
d = 1, 2, 3,4, Hi > 1/2 implies (5.8). When d > 5, then the condition (5.8) is implied by the 
following 

2 7 

dd, > 1 - - - ^ 
d 2 d 

by choosing 70 sufficiently small. Theorem 2 is then proved. □. 
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